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ABSTRACT 



Context. Knowledge of the age distribution of stars in young clusters and associations is crucial to constrain models of star formation. HR 
diagrams of different young clusters and associations suggest the presence of age spreads, but the influence of errors on the derived ages is still 
largely debated. Determination of lithium abundances in low-mass stars represents an alternative and robust way to infer stellar ages. 
Aims. We measured lithium in a sample of low mass members of the young (4-5 Myr) cr Ori cluster with the main goal of investigating its star 
formation history. 

Methods. Using the FLAMES multi-object spectrograph on VLT/UT2, we obtained spectra of 98 candidate cluster members. The spectra were 
used to determine radial velocities, to infer the presence of Ho- emission, and to measure the strength of the Li i 670.8 nm absorption line. 
Results. Using radial velocities, Ha and Li, together with information on X-ray emission, we identified 59 high probability cluster members. 
Three of them show severe Li depletion. The nuclear ages inferred for these highly depleted stars exceed 10-15 Myr; for two of them these 
values are in good agreement with the isochronal age, while for the third star the nuclear age exceeds the isochronal one. 
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1. Introduction 

The time scale of stellar birth within molecular clouds is one 
of the main open issues in star formation theory. Two sce- 
narios have been proposed to explain the observational re- 
sults: Rapid Star Fo rmation (RSF) or Slow Accelerating Sta r 
Formation (SASF) jElmegreenl l200fi IPalla & Stahierj|2002h . 
HR diagrams, which are the main tool to determine Pre- 
Main Sequence (PMS) stellar ages, show a spread of >10 7 
years in individual clusters and associations, supporting the 
SASF model, but the effect of various sources of errors, both 



observational and theoretical, is largely debated (IHartmann 
2001; Bu rningham et al. I l2005h . Lithium (Li) abundances can 
be used as an indepen dent an d robust method to determine 
ages of young objects (IMartfnll 1 998allbl) . since low-mass stars 
in the range 0.5 -0.08 M H deplete the i r initial Li content d uring 
the PMS phase jBodenheimerll 19651; iBaraffeet al]|l998l) . The 
timescale of Li depletion depends on mass, with higher mass 
stars (0.5-0.3 M ) starting to burn it aft er 5-10 Myr and l ower 
mass stars (M<0.2 M Q ) after 20-30 Myr. lPalla etafl §005) em- 



Send offprint requests to: G. G. Sacco, e-mail: sacco@astropa.unipa.it 
* Based on data collected at the ESO Very Large Telescope, Paranal 
Observatory, Chile [program 074.D-0 136(A)] 



ployed the Li age-dating method among low-mass members of 
the Orion Nebula Cluster and found four Li-depleted stars with 
nuclear ages ~ 10 Myr. 



The cr Ori cluster was discovered by ROSAT ( I Walter et al 
1997b around the 09.5 V bin ary st ar cr Ori AB (dis- 



tance 350^ij 6 pc, IPerrvman et al. 19971) . Its low-mass stellar 
population has been intensively studied by photometry and 

low-resolution spectroscopy in the optical and near-infrared 
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(Zapatero Osorio et al. 2002; Barrado v Navascues et al. 2003; 

Sher ry etal.1 12004 ISchoTz& Eisloffell 12004 lOliveira et al j 



2004J), whi le the X-ray properties h ave been investigated most 



recently by Franciosini et al. ( 20061) using XMM-Newton 



Oliveira et al.l (12002b deter mined an isochronal m edian age 



of 4.2+" Myr. Subsequently, lOliveira et all d2004 presented 
evidence for a large age spread (-30 Myr) in the I/I-J Color- 
Magnitude Diagram (CMD), and suggested that part of it 
could be due to p hotometric variability of the PMS stars. 
Kenvon et al. d2005h have found bona fide cluster members 



with small Li line equivalent widths (EW ), possibly indicat- 
ing a certain amount of depletion. Finally. I Jeffries et al.1 (12006) 
have discovered two kinematically separate populations of 
young PMS stars, one concentrated around cr Ori AB sharing a 
common radial velocity with this star (vi = 31.0 + 0.5 km/s), 
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and the second one, more dispersed in the sky, with a radial 
velocity similar to that of the Orion OB la and lb associations 
(v 2 = 23.8 ± 1.1 km/s). 

In this Letter, we report on the discovery of three high- 
probability members of the main cr Ori cluster with Li abun- 
dances a factor of about 1000 below the interstellar value. 
This result was obtained as part of a VLT/FLAMES survey to 
study membership, Li abundances and accretion diagnostics of 
a large sample of K and M stars around cr Ori AB (Sacco et al., 
in preparation). 

2. Observations and Results 

2.1. Target selection and observations 

We have selected 98 cluster candidates from previous stud- 
ies. In order to secure a high fraction of cluster mem- 
bers, we have granted higher pri ority to stars detected in 
X-rays by iFranciosini et al. (2006) and with isochronal ages 
<10 Myr from optical and infrared CMDs. All stars of 
the sam ple have an infrared counterpart in the 2MASS 

catalog (Skru tskie et alj 2006 ). while for 79 of them opti- 
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cal p hotomet ry is available in the literatur e dWolkl 11996; 
Zapatero Osorio et al. 2002; Barrado y Navascues et al. 2003; 



Burningham et al.l l2005l;lKenvon et alf2 005) 



VLT/UT2 dPasquini et al 
MEDUSA mode with 
HR15N grating covering 



The sa mple stars were ob served using FLAMES on 
~* 120021) . FLAMES was operated in 
the high resolution (i/A/l=17000) 
the spectral range 647-679 nm. 
The Field of View (FoV) was centered at RA=05 A 38 m 48*.9 
and Dec=-02 flf 34"' 22 J and is almost coincident with the FoV 
of XMM-Newton (IFranciosini et al. I l2006l) . The observations 
were divided in six runs between October and December 2004, 
for a total exposure time on target of 4.3 hours. Data reduc- 
tion was performed using the GIRAFFE girBLDRS p ipelin e 
1.12, following the standard steps (Bl echa & Simondl 120041) . 
Sky subtraction was performed separately using the average of 
six-seven sky spectra. 

2.2. Membership 

We have measured radial Velocities (RVs) using the IRAFQtask 
FXCOR by Fourier cross-correlation with two template spec- 
tra chosen among sample stars with no accretion signatures. 
For 86 stars we derived a median RV as the weighted average 
of the six values measured for the different runs, while the re- 
maining 12 stars turned out to be candidate binaries and will 
be discussed in the forthcoming paper. Cluster RV was deter- 
mined by fitting the observed RV distribution, shown in Figure 
Q] with the weighted sum of 2 gaussian distributions, one for 
the cluster and the other for the field. The best fit yields a value 
of v c =30.91+0.90 km/s for the cluster and v F =43+ 36 km/s 
for the field. Our cluster velocity is in ex cellent agreem ent with 
that (vc=3 1 .0+0.5 km/s) determined by Jeff ries et al.1 (2006). 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
research in Astronomy, Inc., under contract to the National Science 
Foundation. 
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Fig. 1. Observed RV distribution of all the stars of the sample, 
with overplotted (dashed line) the curve resulting from the fit- 
ting with the sum of two gaussians. 



We considered as cluster members stars with RVs within 
3 cr from the average, yielding 61 members and 25 field stars. 
The statistical contamination of the cluster member sample, 
estimated by integrating the field star distribution between 
Vc - 3cr c and Vc + 3cr c , is ~2 stars. The second requirement 
on membership is the presence of Ha in emission. Among the 
25 RV non-members, 24 are field stars also according to Ho-, 
while one star, which appears young based on Ha, has a ra- 
dial velocity (v=23.73+0.45 k m/s) similar t o that of the sec- 
ond population discovered by Jeffrie s et al. ( 2006b . Note that 
all the non-members have Li I 670.8 nm pseudo Equivalent 
Width (pEW) smaller than 200 mA. Finally, among the 61 RV 
cluster members, 56 are PMS stars also from Ho- emission and 
from the strength of the Li line; two stars with Ha in absorption 
and no Li are most likely field stars. 




Fig. 2. Spectra of the three Li-depleted stars listed in Table 
Q] namely, SWW127 (dotted line), J0539 14.5-022834 (dashed 
line), and SE51 (continuous line). For comparison the spec- 
trum of an undepleted star in same temperature range of those 
reported in Table |2]is also shown (dotted-dashed line) 

The remaining 3 RV cluster members, listed in Table Q] 
have Ha in emission and Li pEW less than 200 mA. As shown 
in Figure |2] the Li line is clearly identified (EW=150 mA) 
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Table 1. Membership indicators of Li depleted stars. 



Object 




Ha EW 


Log L x 


PMS 




(km/s) 


(A) 


(erg/s) 


indicator 


SE51" 


29.5±0.5 


-1.1 


29.22 


variable 


SWW127* 


33.2+0.5 


-2.6 


29.51 




J053914.5-022834 1 ' 


29.2+0.5 


-4.6 


30.17 


NIR ex. 



Zapatero Osorio et al. (2002). 



in the spectrum of SE51, but not in those of SWW127 and 
J0539 14.5-022834. These stars have an X-ray counterpart and 
result to be PMS stars in CMDs. Moreover J0539 14.5-022834 
shows excess in K-L dOliveira et alJl2004 and SE51 is char- 
acterized by photometric variability ( Scholz & E isloffel 2004). 
Considering that only two out of 22 field stars in our sample 
have an X-ray counterpart (2 field stars are out of the XMM 
FoV) and that the probability of having a field star with RV 
between vc - 3<xc and vc + 3<x c is 0.056, the probability of 
finding, among the RV cluster members, 3 field stars with an 
X-ray counterpart is less than 5xl0~ 3 . We conclude that, al- 
though physical membership must be confirmed by, e.g., proper 
motion studies, we regard as unlikely that all of them are non 
members. 



els of Baraffe et al.l d 1998b . The resulting Li abundances are 
listed in Tableland are a factor of > 1000 below the interstel- 
lar value (A(Li)=3.3). We stress that even relatively large errors 
in T e ff and/or log g would not greatly affect Li abundances (or 
depletion factors), given the extremely small pEWs (or upper 
limits to pEWs). Similarly, the low abundances cannot be ex- 
plained by veiling, since it would imply a value of r » 1, or 
by other sources of error, such as unce rtainty in the pseudo - 
continuum. Contrary to the suggestion of lKenvon et al.l (12005b . 
we believe that the measured low pEWs represent real Li de- 
pletion, since Pavlenko (2001) models provide a correct treat- 
ment of the Li line formation in cool objects and pEWs below 
~170 mA imply that only TiO contributes to the 670.8 feature. 
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2.3. Lithium depleted stars 

pEWs of the Li line were measured using the IRAF task 
SPLOT, by direct integration of the line profile over an interval 
of 0.2 nm. Measured pEWs may be affected by spectral veiling. 
We have estimated r, the ratio of the excess emission to the pho- 
tospheric continuum, for all stars of the sample from the mea- 
surement of the EWs of three absorption lines (V I 662.5 nm, 
Ni I 664.3 nm, Fe I 666.3 nm) in their spectra and in those 
of 1 1 unve i led co mparison stars with similar temperatures (see 



Palla et al.l (120051) ). In those cases where we were not able to 
measure at least 2 lines, because of low S/N or high veiling, we 
have considered r indeterminate and the measured Li pEWs are 
lower limits to the true values. Derived r values range between 
and 1.4, with errors of about 0.1-0.2. Figure [3] shows pEWs 
of the Li line, corrected for veiling, for stars with available op- 
tical photometry. The median pEW is 590 mA, with a typical 
dispersion of ~100 mA. The scatter in Li pEWs could be due 
to measurement errors, but we cannot exclude the presence of 
some partially depleted stars. The three stars reported in Table 
[TJare not veiled (r=0.0-0.2) and their Li pEWs are less than 200 
mA, indicating a large amount of Li depletion. 

We have determined Li abundances for these thre e stars us- 
ing Cu rves of Growths of Pa lla et al. (2006), based on Pavlenkol 
d200ll) models and spectral code. Surface gravity was fixed at 
log g=4. dex, the expe cted value for a 4 Myr old late type star 
by iBaraffe et al. ([l998) models. The effective temperature of 
the depleted stars is given in Table|2] in the case of J053914.5- 
022834 it has been derived fr om the spectral type measured by 
Zapatero Osorio et al.ld2002h . while T eff of SE51 and SWW127 



are estimated from the color indexes R-I and I-J, using the mod- 



es 1.0 1.2 1.4 1.6 1.8 2.0 

R-I 



Fig. 3. pEWs of the Li line corrected for veiling as a function of 
R-I color for cluster members with available optical photom- 
etry. Upward and downward arrows denote lower and upper 
limits. 



3. Discussion and Conclusion 

The three Li depleted, high probability members of the cr Ori 
cluster represent an example o f post T Tauri stars (following 
the definition of Martin lll997h i. and allow us to investigate its 
star formation history from a novel point of view. Isochronal 
and nuclear ages (Ihrd) and masses (Mhrd) for the three Li- 
depleted stars are listed Table [2] These values have been de- 
rived from the models of Palla & Stahler ll 19991) using the T eff 
and luminosity range given in the table. In the last two columns 
we list the nuc l ear m ass (M Li ) and age (tjj) estimate d follow- 
ing IPalla et all d2005l) . based on lBildsten etail dl997l) . The lat- 
ter provides analytical prescriptions to derive the age vs. stellar 
luminosity and the age vs. Li abundance relations for a fully 
convective star undergoing gravitational contraction at a nearly 
constant T e ff and assuming fast and complete mixing. In Figure 
|4]we display the mass vs. age relations at fixed luminosity (pos- 
itive slope) and Li abundance (negative slope) predicted by the 
models. The intersection of the two curves gives the mass and 
nuclear age of the star. 

We find that the age and mass of SE51 and SWW127 de- 
rived from Li are in very good agreement with the isochronal 
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Table 2. Li depleted stars and their properties. 



Object 


Teff (K) 


Log(L/L Q ) 




LiEW 


A(Li) 


M H rd 


tHRD 


M Li 


tu 




(RJ/IJ) 


(RJ/IJ) 




(mA) 




(Mo) 


(Myr) 


(Ma) 


(Myr) 


SE51 


3750+20/3750+60 


-0.91+0.01/-0.9L 


r0.02 


150+50 


-0.5+0.5 


0.52-0.53 


12-13 


0.49 + 0.04 


10.5 ±0.6 


SWW127 


3645±15/3720±60 


-0.92±0.01/-0.94=i 


=0.02 


<90 


<-1.0 


0.42 - 0.49 


8.4-14.2 


> 0.44 


> 10 


J053914.5-022834 


3345+70 


-L14±0.06 




<100 


<-0.7 


0.23 ± 0.05 


4.5 ± 1.5 


> 0.4 


> 15 



values obtained using both R-I and I-J colors. On the contrary, 
the measured amount of Li depletion for J0539 14.5-022834 is 
too large for the isochronal age and low mass from evolutionary 
tracks. 



at high resolution to measure lithium abundances on a larger 
sample of stars can help to fully characterize the star formation 
history of the cluster. 
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Fig. 4. Mass vs. age plot for the Li depleted stars SE51 
and SWW127. T e ff and bolometric correction have been de- 
rived using the color indices R-I (left- hand panels) and I- J 
(right-hand panels) and the models of iBaraffe et alj dl998l) . 
Continuous curves represent the locus of constant luminosity 
(positive slope) and constant Li abundance (negative slope), 
while dashed curves indicate the uncertainty range. The circles 
with the error bars give the mas ses and ages computed from 
models of Pall a & Stabler! d 1999b . 



Considering that most of the remaining stars of our sam- 
ple have pEWs consistent with the interstellar value, our re- 
sults support the view that the bulk of the <x Ori population has 
an age < 4-6 Myr. However, three low-mass, high probability 
members are definitively older than ~10 Myr. Therefore, we 
propose that cr Ori has been forming stars on a time scale of 
> 10-15 Myr, as found in the case of the Orion Nebula Cluster. 

Finally, we note that our sample has been selected requiring 
an isochronal age from CMDs less than ~10 Myr. This bias has 
precluded the identification of additional Li-depleted, old stars 
that might exist in the cr Ori cluster. Further observations at 
low spectral resolution to derive stellar parameters, as well as 
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